Ab initio estimate of temperature dependence of electrical conductivity in a model 
amorphous material: hydrogenated amorphous silicon 
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We present an ab initio calculation of the DC conductivity of amorphous silicon and hydrogenated 
amorphous silicon. The Kubo-Greenwood formula is used to obtain the DC conductivity, by thermal 
averaging over extended dynamical simulation. Its application to disordered solids is discussed. The 
conductivity is computed for a wide range of temperatures and doping is explored in a naive way 
by shifting the Fermi level. We observed the Meyer-Neldel rule for the electrical conductivity with 
£'mnr=0.06 eV and a temperature coefficient of resistance close to experiment for a-Si:H. In general, 
experimental trends are reproduced by these calculations, and this suggests the possible utility of 
the approach for modeling carrier transport in other disordered systems. 

PACS numbers: 72.80.Ng, 72.20.-1,71. 23.Cq 



I. INTRODUCTION 

Amorphous semiconductors are among the most im- 
portant electronic materials. a-Si:H is the universal 
choice for TFT applications in laptop displays and im- 
portant photovoltaic material. Thin films of a-Si:H 
are one of the promising active elements for uncooled 
microbolometersi as employed in focal plane arrays for 
night vision applications. Here, the temperature (T) de- 
pendence of the electrical conductivity, and the Temper- 
ature Coefficient of Resistance (TCR) is of paramount 
importance. Experiments suggest that there are well- 
defined conductivity regimes: variable-range hopping at 
the lowest temperatures, phonon-induced delocalization 
at higher temperatures and ultimately a metallic form 
of conduction (albeit with strong scattering). The ven- 
erable Meyer-Neldel rule^ii^ (MNR) correlating the expo- 
nential prefactor of the DC conductivity with the activa- 
tion energy has been observed to be almost universally 
valid in disordered systems. A microscopic understand- 
ing of these effects is highly incomplete, though signifi- 
cant progress has been made on MNR'^'^. 

Disordered solids pose definite challenges for transport 
modeling. The Boltzmann equation is best suited to long 
mean free paths, and is difficult to apply in these sys- 
tems with localized states and defects. Moreover, the 
topological disorder must have a significant influence on 
the conduction and a realistic calculation of transport 
must include this disorder explicitly (eg through a cred- 
ible structural model). 

Modern density functional simulations of materials 
routinely provide the Kohn-Sham"" eigenvalues and or- 
bitals at each step in a thermal Molecular Dynamics 
(MD) simulation. The natural approximate connection 
of these quantities to electrical transport is provided by 
the Kubo-Greenwood formula (KGF)^-^. The first ap- 
plication of the KGF to an atomistic model with com- 
puted electronic structure was by Allen and Broughton^ 
for (metallic) liquid Si with a tight-binding Hamiltonian. 
In this paper, we explore the utility of the KGF with 
Kohn-Sham spectral properties to estimate the tempera- 



ture dependence of the conductivity in a-Si and a-Si:H. 
An adiabatic approximation is used for modeling the in- 
elastic scattering of the phonons (by thermally averaging 
the KGF over a constant T MD trajectory). The under- 
lying idea is that if the thermal fluctuations drive the un- 
occupied and occupied states close together (in the spirit 
of Landau-Zener^ tunneling), then the energy conserving 
6 functions are nonzero, and a finite contribution to the 
conductivity is made by such instantaneous "snapshots" 
if the dipole matrix element is non-zero. Doping, the 
MNR and TCR are obtained and shown to be in reason- 
able agreement with measurements in a-Si:H. The results 
are sufficiently encouraging to justify a fuller exploration 
of the applicability of the method to disordered systems 
more generally, including other amorphous semiconduc- 
tors and glasses and conducting polymers. 

The paper is organized as follows. In Section II, we 
discuss the underlying logic of the approach and some 
salient previous calculations. In Section III, we briefly 
describe models and simulation procedures employed. In 
Section IV, we present results for a-Si and a-Si:H includ- 
ing the T-dependent conductivity, the effect of doping on 
the conductivity and a discussion of MNR and TCR for 
a-Si:H. Finally in Section V we draw conclusions. 



II. APPROACH 

A. Overview 

Previous work on amorphous semiconductors and 
glasses has shown that Kohn-Sham eigenvalues conju- 
gate to localized Kohn-Sham states are specially sus- 
ceptible to the motion of the latticeiSiiii (thermal fluc- 
tuation in the atomic coordinates involved in localized 
states leads to a strong modulation of such eigenvalues 
and eigenvectors). Indeed, the RMS fluctuation of eigen- 
values conjugate to localized eigenvectors is proportional 
to the localization of the eigenvector as measured by in- 
verse participation ratio^i. Localized states appear in 
the gap, or in the spectral band tails. Since they are 
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close to the Fermi level, they play a role in transport. 
The most elementary link between these states and ener- 
gies and the conductivity is the KGF, which is a natural 
approach to strong scattering systems with short mean 
free paths. With this in mind, we use the Kubo formal- 
ism and current ab initio methods in conjunction with 
fairly realistic (structurally plausible) supercell models of 
a-Si:H to probe the temperature dependence of the DC 
conductivity. The novelty of our work lies in 1) using 
Kohn-Sham states and energies in Kubo's expression for 
the conductivity and 2) adopting a Born-OppenheimeriS, 
like approximation of thermally averaging the KGF (us- 
ing instantaneous atomic configurations obtained in the 
course of a thermal simulation at constant temperature 
as "snapshots"). Such an approach certainly has limita- 
tions: for low temperatures in which one expects variable 
range hopping between defects the correct link to first 
principle simulation would seem to require solutions of 
the time-dependent Kohn-Sham equations^'*'^^ or an ap- 
proach more akin to a Miller- Abrahams^ - model of the 
conductivity. There are several studies which use the 
KGF to compute the static lattice conductivity of amor- 
phous materialsiSiiiii^. The computed conductivity van- 
ishes for localized states in a static lattice. 

As an alternative to the weak scattering theory of 
Ziman"^^ for liquid metals, investigators have used the 
KGF with variants of Car-Parrinellc^ simulation to ex- 
plore the conductivity of liquid metals including for ex- 
ample, 1-NaSn alloys^" and 1-Na^^, with quite satisfac- 
tory results. It was found that Brillouin zone sampling 
had to be treated with care, which is reasonable for a 
metallic system. The success of such calculations is one 
motivation to extend the approach to doped amorphous 
semiconductors^. Another significant development of re- 
cent years is the Keldysh Green's function approach of- 
fered to properly handle the effects of contacts and finite 
potential drops for molecular electronics. An example of 
this approach is the code TRANSIESTA^a. 



B. Connection to many-body formulation 

The KGF has been applied on many occasions, and in 
liquids in a mode very similar to what we report here. 
However, there is usually no discussion of the underlying 
assumptions that connect the original many-body formu- 
lation of the KGF to its usual use with MD averages and 
single particle states. For example, the KGF as used 
in these computations is not formulated for use with in- 
elastic processes, and thus the use of MD trajectories and 
averaging requires some discussion, which we offer in this 
section. 

Consider a periodic perturbation Hint 

H„,t ^ Fe-'"^^ + F*e"^* (1) 

acting on a system described via a many-body (electron 
and phonon) Hamiltonian H. The transition probability 



from initial state '^i to final state 4"/ in the interval (Ivf 
under the action of an external field (Hint) is^^ 

— \{^^Jf\F\^^J,)\^SiEf-E,-nu;)duf (2) 

where ^'^ are eigenstates of many-body Hamiltonian 
H 

H-^f^Ef^f, H ^ Hq + H'^-P'' . (3) 

Here, Ef and f may be estimated by the ordinary time- 
independent perturbation theory from the eigenfunctions 
^fp"' and eigenvalues Ej^^ of 

Ho = He + Hion ■ (4) 

We adopt a Born-Oppcnheimer description 

iJo4.(o)^ij(o)^(o) (5) 

with, 

*f = <i>e,e„^ (6) 

and <i>ep is a many electron wave function, Qy^ is a many 
phonon wave function. The perturbation solution of Q 

is 

*, = ^a.,vl/(o), a,, = S,, + a^J+a^J + --- (7) 

E.^E^^+E^^+eP + ... (8) 

Detailed expressions can be found in standard textbooks. 
The coupling H^~p'^ between electrons and phonons need 
not be small, in principle eigenfunctions and eigenvalues 
of H can be calculated to any order of H'^~p^ . 

In — > limit, we only keep the interaction F^ be- 
tween electrons and external field. Eq. [5] is modified to: 

GpGq Up 

X S[{EfUEf+Ef^ + ---) 

-(i?f)+i?P+i?P+...)-Md^/ (9) 

Eq. [HI shows that for a; = 0, the conservation of energy is 
for the whole (electron-|-phonon) system. Phonons may 
assist the transition between two many-electron states 
with different energy if the phonons can serve as source 
or sink of energy. 

At finite temperature, initial state can be any one 
of all possible states. To calculate the total absorption 
power in a sample, we average initial state with Boltz- 
mann weighting, and sum over all possible final states. 

The Born-Oppenheimer approximation (BOA) starts 
with a stable lattice Aq : {X°, Y°, , ■ ■ ■ X%,Y^, Z%) as 
its zero order configuration. The sum over all possible 
phonon states Vp, one phonon, two-phonon, •■•etc., in 
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Eq. [9] means that we explore all possible configurations 
of the lattice around our original conformation Aq. In 
the configuration space spanned by all N nuclear coordi- 
nates {Xi,Yi, Zi, ■ ■ ■ XpfjYpf , Zn), those configurations 
included in in Eq. [9] form a 3N dimensional region 
V around Aq. 

The next step in the transition to current thermal sim- 
ulations is the treatment of the ions as particles with tra- 
jectories governed by the classical equations of motion. 
For a system like a-Si:H, this is a defensible approxima- 
tion, at least above the Debye temperature. Now, con- 
sider a MD simulation commencing from a configuration 
in the neighborhood of Aq. One would then approximate 
Eq. [Hby 

(10) 

where n is the number of MD steps. The assumption here 
is that the MD trajectory faithfully reproduces the dy- 
namical processes implicit in Eq. [2l In each MD step, the 
Kohn-Sham scheme gives the fully dressed single-particle 
energy for the electron-phonon coupling of the last MD 
step. In a sense, each MD step can be interpreted as a 
sum of one specific series of Feynman graphs in set P, if 
we classify Feynman graphs according to the configura- 
tion space points in 2?. 

For a chosen basis set, one can expand many-electron 
wave function $ as a linear combination of Slater deter- 
minants of single particle functions "tpr, 

= e,i,....,„4^ct^---ct„|000---0). (11) 

qiq2---qN 

The interaction between electrons and external field E 
along the x direction can be written in single particle 
form as: 

Fe ^ {^)Y.mx\^s)ctcs (12) 

rs 

from which the usual Kubo-Greenwood formula is ob- 
tained as we discuss in the next subsection. 

Depending upon the order of perturbation of Eq. 
[71 Eq. [5] describes various many-phonon assisted pro- 
cesses. This agrees with the Green function formulation. 
Green function metho d^^i^^ starts from free single elec- 
tron propagator and free single phonon propagator, this 
is contrast with Kubo formula which start with exact 
many-body wave function for the whole electron-|-phonon 
system. Technically, the former will be much easier than 
the later. Both acoustic and optical phonon assisted 
hoppingSS are worked out and applied to the transition 
rate from one localized state to another localized state, 
and the Meyer-Neldel relation is viewed as a consequence 
of assisted activation whenever large activation energies 
compared to typical excitation are involved^iSl. 



C. Application of KGF to disordered solids at 
finite temperature 

The derivation of the KG formula from linear re- 
sponse theory and the fluctuation-dissipation theorem 
is available in the original literature^, and elementary 
derivations are provided in standard books on transport 
in amorphous systems^i^. In these latter derivations, 
first-order time-dependent perturbation theory (Fermi's 
Golden Rule) is employed to deduce an expression for the 
AC conductivity cr{uj). From either derivation, one ex- 
pects the KGF to be valid in the weak-field limit and for 
elastic scattering processes. The form for the diagonal 
elements of the conductivity tensor for the static lattice 
that emerges from the preceding discussion is: 

ni 

X (5(e„ - Si - huj) (13) 

where a = x,y, z, fj;- is the Fermi distribution, e and m 
are the electronic charge and mass, pa is a component 
of the momentum operator, ijii and £i are the eigenstates 
and eigenvalues; Q, the cell volume. The AC conduc- 
tivity is then (t(w) = ^ 'Ylia '^aa- I'^ the rest of this pa- 
per, single-particle Kohn-Sham states and eigenvalues are 
used for the {V'}, {e}- 

In the DC-limit (uj ^0 ) the conductivity takes the 
form 

a(T) = -IJ2 r <yUe)^^de (14) 

a ^ 

where: 

ni 

(15) 

We include thermally-induced electron state and energy 
fluctuations near the gap, by averaging expressions such 
as the preceding over a thermal simulation. The DC con- 
ductivity is computed from trajectory averaged quanti- 
ties such as: 

^c.a(e) = ^ E MnWm?5{ei-e)5{e\-e), 

ni 

(16) 

where the bar denotes the average and we emphasize 
the dependence of the various terms on the simulation 
time t. This average then picks up thermal broadening 
effects in the density of states, and also include time- 
dependence in the dipolc matrix element. A Gaussian 
approximant for the 5 function with a width of 0.05 eV 
is used in our calculations. We have repeated our calcu- 
lations for widths between [0.01-0.1] eV, and our results 
do not change appreciably. We insert thermally averaged 
diagonal elements as in Eqn. [11] into Eqn [H] to obtain 
the conductivities reported in this paper. 
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The KGF includes three possible transitions: (i) lo- 
calized state to localized state; (ii) localized state to 
extended state; (iii) extended state to extended state. 
For amorphous semiconductors, the upper tail of valence 
band and the lower tail of conduction band are localized. 
Below room temperature, carriers are mainly distributed 
in localized states, and (i) is dominant conduction mecha- 
nism. With increasing temperature, there are more holes 
in the upper valence tail, which are available for the elec- 
trons in lower side of mobility edge to transit. In addi- 
tion, the activation energy is lowered by thermal fluctua- 
tions. At higher temperature, (ii) becomes more impor- 
tant until there are enough carriers in extended states, 
and (iii) will becomes the most effective conduction chan- 
nel. 

It is believed?^ that the temperature range from 
ca 300K-700K is dominated by "phonon-induced 
delocalization"— , in which the phonons aid the hop- 
ping (though not to the extent of polaron formation). 
Phonon-induced delocalization has been computed from 
the time-dependent Kohn-Sham equationi^. It was ob- 
served that if electron energies became close (because of 
phonon-induced variation in energies), mixing between 
the nearly degenerate states could occur, which results 
in delocalizationi^. This leads to a monotonic and irre- 
versible decrease in electron localization. In the present 
calculation, the localization fluctuates (see Fig. 2 of Ref. 
[Tsh . Thus, our calculation ignores the explicit phonon- 
driven delocalization except in so far as this effect is cap- 
tured in the thermal averaging scheme. 

There are sources of error in any calculation of this 
type, which we must mention: (1) We use the Kohn- 
Sham states computed in the LDA as input into the 
KGF. (2) We employ very small structural models (ca 
70 atoms); (3) Except for a few test cases, we have 
used a minimal (singie-zeta) basis. In fact we do this 
in part because this approximation gives a gap close to 
the physical gap in a-Si\ (4) Our simulations run for fi- 
nite time, though we have taken some care to verify that 
the reported conductivities are adequately converged; (5) 
There are physical mechanisms that should broaden the 
density of states. Two of these arise from Brillouin zone 
integration and the sparse and discrete sampling of the 
tail states due to the small cell size. (6) Classical dy- 
namics (no quantization of the vibrations) is employed 
throughout. These errors are largely systematic, and as 
such, difficult to measure. 

On the other hand, there are conditions that amelio- 
rate the situation. Where point (1) is concerned: quasi- 
particle calculations for crystalline Si computed in the 
GW approximation and the Kohn-Sham states are nearly 
identical"^^-. Also, there is a strong thermal averaging ef- 
fect: even at room temperature, localized tail states can 
fluctuate by tenths of an eV (far greater than fcT); this 
helps significantly in connection with points (2), (3) and 
(5) above. Finally, there is a convincing body of empir- 
ical evidence that the KGF employed as we do here (or 
with more primitive Hamiltonians) provides reasonably 



quantitative predictions for the DC conductivity. The 
use of a SZ basis is convenient computationally, and has 
the helpful feature that the gap is much closer to exper- 
iment than a complete basis calculation. 



III. HAMILTONIAN AND STRUCTURAL 
MODELS 

The ah initio local orbital code SIESTA^- was used to 
perform the density functional calculations. We used a 
local density approximation for the exchange (LDA) us- 
ing the Perdew and Zungcr expression'^^ . Norm conserv- 
ing TrouUier-Martins pseudopotentials'^^ factorized in the 
Kleinman-Bylander form^^ were used. We used a min- 
imal single C basis set for both Si and H^.We solved 
the self consistent Kohn-Sham equations by direct diag- 
onalization of the Hamiltonian and a conventional mixing 
scheme. We used the T{k ~ 0) point to sample the Bril- 
louin zone in all calculations. No scissors correction was 
used since the SZ optical gap is close to experiment (a 
scissors shift is necessary for a polarized basis). 

We have used two different well relaxed models of aSi^ 
(64 Si atoms >22 and aSieiHio (61 Si and 10 H atoms )^ 
models. We have prepared these models for six differ- 
ent temperatures 200K, 300K, 500K, TOOK, lOOOK, and 
1500K. In each cases we followed the following proce- 
dures. The two models were annealed to a particular 
temperature for 1.5 ps which is followed by equilibration 
for another 1.5 ps. Once the models are well equilibrated, 
we performed a constant temperature MD simulation for 
another 500 steps to obtain an average DC conductivity 
for the respected models at a given temperature. 



IV. RESULTS 
A. Amorphous silicon: a-Si 

1. Conductivity 

We have studied the electrical properties at different 
temperature by using the inverse participation ratio (I) 
which is defined as 

N 

where N is the total number of atoms and qi{s) is the 
MuUiken charge residing at an atomic site i for an eigen- 
state with eigenvalue e with X^i^il^il^)] ~ ^- I unity 
for an ideally localized state and 1/A^ for an extended 
state. 

In Fig. [l]we show the instantaneous IPR as a function 
of energy for six different temperatures. As the tem- 
perature increases from 200K to 1500K, the optical gap 
is reduced and eventually at higher temperature all the 



states become extended with no energy gap in the den- 
sity of states. The gap closes both because of thermal 
fluctuations on the eigenvalues and as a harbinger of the 
transition to a metallic state at sufficiently high temper- 
atures. We emphasize that these are only snapshots: the 
instantaneous IPR of a well-localized state can varj*!^ by 
a factor of « 2. This is understood as a consequence 
of the fluctuation of the eigenvalues fas eigenvalues ap- 
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FIG. 1: Instantaneous snapshot of IPR (I) versus energy for 
different temperatures in aSi^i model. The Fermi level falls 
in the gap near E = —3.8eV. Localization is significant near 
the band edges, and note that the gap closes at higher tem- 
perature: at 1500K the cell has melted, at lOOOK the gap has 
closed because of thermal fluctuations in the eigenvalues. 
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FIG. 2: (color online) DC conductivity of intrinsic aSie4 av- 
eraged over 500 configurations computed at different temper- 
atures. The solid symbols are from our work, open symbols 
are from experiment: open square from Ref. |39| . open dia- 
mond from Ref. I8i , and open triangle from Ref. |4l|. The 
error bars are from equilibration time and discrete sampling 
of the density of states. 



where Ea is the activation energy {Ea = Ec — Ep or 
Ea = Ep — Ey) and CTq is the pre-exponential factor 
of the conductivity. By dividing the dc conductivity in 
two regions of low temperature (T < A50K) and high 
temperature (T > A50K) we extracted the Ea and Co- 
For low T, we have obtained Ea ~ 0.34 eV and CTo ^ 
4 f]-icm-i. For high T, Ea 0.45 eV and do - 1 x 
10^ n-^cm^\ 



The DC conductivity of aSi^i is accumulated over 500 
instantaneous configurations for different temperatures: 
T=200K, 300K , 500K, TOOK, lOOOK, 1500K, and 1800K. 
At a temperature of 1800K, the system is actually a liquid 
with a diffusion coefficient of D ^ 1.6 x 10"'* cm^s~^ 
and a DC conductivity of ~ 0.3 x 10''' ^l~^cm~^ which is 
reasonably close to the measured value of (1.0 — 1.3) x 
10* n-^cm-^M and value of 1.75 x 10'' n-^cm-^ obtained 
in another computation.'"'. 

In Fig. [2] we have shown the DC conductivity of aSigi 
as a function of temperature. The results from experi- 
ment for selected temperatures are also shown. As we 
can see from Fig. [1] increase in the temperature of the 
system enhances delocalization of the states and eventu- 
ally closes the optical gap to change the material from 
semiconductor to metal. In doing so the DC conductiv- 
ity changes from 0.31 x lO"*" f^-^cm-^ for T=200K to 
0.24 x 10^ n-^cm-^ for T= lOOOK. 

The temperature dependence of DC conductivity can 
be written as 



B. Doping 

It is known that doping and temperature change result 
in a shift in the position of Fermi level within the optical 
gapi^^. In our simulation, we have computed the DC 
conductivity for a given doping by shifting the Fermi level 
from its intrinsic position towards the conduction band 
edge or valence band edge in steps of O.leV. This pro- 
cedure allows us to "scan" the optical gap and compute 
conductivity for different doping levels of n-type as well 
as p-type. Our scheme is quite different from experiment. 
For example, it is known that the energy of defect states 
depends upon the location of the Fermi level'''', so that 
a proper calculation of doping (meaning with explicit 
substitution of the donor or acceptor species) should be 
carried out self-consistently. Our procedure is a highly 
idealized version of the "doping problem" , which shows 
how the conductivity varies for without compensation ef- 
fects. Because of these effects, it is not easy to correlate 
the conductivities we predict for a specific Fermi level 
position with the experimental concentration of dopant 
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atoms. We are presently extending this work by explicitly 
including B and P impurities in the system and comput- 
ing the conductivities separately for each doped model 
system.— 

The computed DC conductivity for different temper- 
atures as a function of chemical potential is shown in 
□ 
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FIG. 3: (color online) DC conductivity of aSun as a func- 
tion of doping, averaged over 500 configurations computed at 
difi^erent temperatures versus chemical potential. 



C. Hydrogenated amorphous silicon: a-Si:H 
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1. Conductivity 

In the same way we that analyzed aSi^Ji in the pre- 
vious section, we have started our analysis of aSieiHio 
by computing its electronic properties from the inverse 
participation ratio. In Fig. 21 we have shown the IPR I 
of aSieiHiQ for different temperatures. As can be seen 
from the figure, the optical gap decreases with increasing 
the temperature which is attributed to phonon induced 
delocalization and structural rearrangements. 

The DC conductivity of aSieiHiQ as a function of tem- 
perature is shown in Fig. [5] with comparison from ex- 
perimental results from Beyer et aZ4^. As we can see 
from Fig. \^ increase in the temperature of the system 
enhances delocalization of the states and eventually elim- 
inating the optical gap to change the property of the ma- 
terial from semi-conductor to metal. In this case, the 
DC conductivity changes from 0.54 x 10~^° r2~^cm~^ for 
T=200K to 0.83 x 10^ H-^cir-^ for T= lOOOK. 



FIG. 5: (color online) DC conductivity of aSieiHio averaged 
over 500 configurations computed at difi^erent temperatures. 
The solid symbols are from our work, open symbols are from 
experiment Ref. 46,. The error bars are from equilibration 
time and discrete sampling of the density of states. 



For low T, we have obtained Ea ^ 0.31 eV and 
CTo - 5 n-^cm-\ For high T, Ea 0.36 eV and 
(To 5 X 10'^ ri^^cm^^. These results are in a reasonable 
agreement with the experimental results of Kakalios et 
al^. In studying doped a-Si:H, Kakalios et al. showed 
that for low T, the Ea ranges from 0.16 to 0.21 eV with 
(To - (5 - 10) 17-icm-i. 

The computed DC conductivity for different temper- 
atures as a function of chemical potential is shown in 
Fig. [HI As in the case of aSiQ^, the DC conductivity 
increases as the Fermi energy shifts toward either the 
valence or conduction band from mid gap. At higher 
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FIG. 6: (color online) DC conductivity of aSieiHio as a func- 
tion of doping, averaged over 500 configurations computed at 
difi'erent temperatures versus chemical potential. 

As the temperature increases from 300K to lOOOK, the 
contribution from the matrix elements decreases while 
the contribution from the density of states near the Fermi 
level increases. The temperature dependence of the DC- 
conductivity is most affected by the density of states 
which is considerably broadened for higher temperature. 

2. Meyer-Neldel rule 



in which case the number of available paths or "chan- 
nels" is T-dependent and certainly increasing, reflecting 
an cntropic increase as discussed in MEE^*'. 

For pre-exponential factor (Tq and activation energy 
Ea, the MNR may be expressed as: 

ao = CTooe^''/^MNR (19) 

By performing a linear fit on the DC conductivity re- 
sults, we identified the intercept at (1/T) = to CTq and 
the slope to the activation energy Ea ■ There are number 
of experimental results on a-Si:H which show this expo- 
nential behavior with i?]y[j,j]^=0.067 eV— By plotting 
tj as a function of 1/T for various dopants (n-type as 
well as p-type) we extracted ctq and Ea for aSieiHio and 
the results are shown in Fig. [71 Our calculation gives an 
exponential behavior of CTq as a function of Ea, reflecting 
the Meyer-Neldel rule, with £^mNR= 0-060 eV. 




In a band picture of conduction, one expects an expo- 
nential conductivity a = aoexp{— Ea / kT) . In disordered 
systems (not limited to a-Si:H) there is a more compli- 
cated T-dependence, and in particular a "kink" in the 
conductivity in the vicinity of 400-500K— . The first in- 
terpretation for this kink was the existence of distinct 
low-T and high-T conduction mechanisms with signifi- 
cantly different activation energies (slopes)^. It is ar- 
gued that this interpretation is unlikely, and that the ef- 
fect arises from the temperature dependence of the Fermi 

48 

energy^. 

For activated processes with activation energy greatly 
exceeding the characteristic excitation energies available 
to the system, it is clear that a fluctuation involving many 
small excitations will be needed to push the system over 
the barrier. The larger the number of ways the neces- 
sary fluctuation may be obtained the more likely it is for 
the process to occur. This line of thinking led to the 
concept of "Multiple Excitation Entropy" (MEE), which 
has clarified the origin of the Meyer-Neldel Rule (MNR) 
in a great many different activated processes^. One can 
think of the KGF as representing a sum over pathways. 



FIG. 7: (color online) The pre-exponential factor as a function 
of activation energy (Meyer-Neldel rule) for aSieiHio. The 
dashed line is an exponential fit which represent Eq. (|19|) . It 
is interesting to compare this result of simulation to a similar 
plot for a large number of experimental samples as in Fig. 7.3 
(Reference [291) . 

The preceding shows that the mechanism(s) respon- 
sible for MNR are present in our simulations. The en- 
hancement of the conductivity with increased activation 
energy may be qualitatively understood in our picture 
as being due to the increase in electron-lattice coupling 
with increasing activation energy (and therefore local- 
ization). Since localized states possess an "amplified" 
electron-phonon couplingiSiii, some compensation is to 
be expected. In our picture, the MNR arises because 
the phonons treat electrons with different localization (or 
Ea) differently, and the effect runs in a direction consis- 
tent with experiment: for doping into the more localized 
states, the electron-phonon coupling is larger and serves 
to modulate the energies more strongly than for more 
weakly localized states with smaller Ea- 
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3. Temperature Coefficient of Resistance 

The other fundamental characteristic of aSi:H is its 
high temperature coefficient of resistance, which makes 
it a candidate for uncooled microbolometer apphcations. 
The temperature coefficient of resistance (TCR) is de- 
fined as 

TCR=-£;^ (20) 
PoT -To 

where p is a resistivity at any given temperature T and 
Po is a resistivity at a reference temperature To (usuaUy 
room temperature). The computed result of TCR with 
a definition of Eq. [201 using To = 300K for aS'igi-ffio 
is shown in Fig. [5] The experimental TCR near room 
temperature is —2.7%K~^ for a-Si:H—. Our calculations 

predict a TCR value of 2.0%K'^ at T = 350^ which 

is in agreement with the experiment. Close to Tq the 
value of TCR is very sensitive to temperature and has a 
wide range of values —(2.0 — 5.0)%K~^. 








Temperature (K) 



FIG. 8: The temperature coefficient of resistance (TCR) for 
aSieiHio as a function of temperature. 



V. CONCLUSION 

We have presented a study of transport in an amor- 
phous material. We used Kubo-Greenwood formula for 
computing the DC conductivity of aSiei and aSigiHio 
for different temperatures. We have also presented the ef- 
fect of doping on the DC conductivity. As the Fermi level 
approaches either the conduction edge or valence edge we 
observe an increase in the DC conductivity. Once the Ef 
exceeds the "mobility edge" we observe a weak tempera- 
ture dependence on the DC conductivity. Though it re- 
quires further investigation (by using various dopants in 
aSisi and aSieiHiQ ), we observe the Meyer- Neldel rule 
with exponential behavior for the pre-exponential factor 
(Tq. The computed result for TCR is in good agreement 
with the experiment. Further study of this method in- 
volves using a richer basis set, and fuller k-point sampling 
in the Brillouin zone (or essentially equivalently, larger 
cells). Further "deconstruction" of the KGF will be un- 
dertaken to explore the instantaneous configurations that 
provide significant contributions to a. 

Our work shows that the simplest implementation of 
the KGF shows promise for computing the electrical con- 
ductivity of disordered semiconductors, and this success 
hints that it may be used with profit on other systems 
beside. 
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